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Tunneling measurements have been carried out on heavily underdoped and slightly overdoped 
Bi2Sr2CaCu208+a: (Bi2212) single crystals by using a break-junction technique. We find that in- 
plane tunneling spectra below Tc are the combination of incoherent part from the pseudogap and 
coherent quasiparticle peaks. There is a correlation between the magnitude of the pseudogap and 
the magnitude of the superconducting gap in Bi2212. We find that the quasiparticle conductance 
peaks are caused by condensed solitonlike excitations. 
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Recent intrinsic c-axis tunneling data obtained in 
Bi2Sr2CaCu208+x (Bi2212) mesas show that the pseu- 
dogap (PG) is a normal-state gap, and the PG and the 
superconducting gap (SG) coexist below Tc Thus, 
the PG in Bi2212 arises either from charge-density waves 
(CDW) or from local antiferromagnetic (AF) correlations 
[or spin-density waves (SDW)] [|l],||j|]. There is a consen- 
sus on doping dependence of the PG in cuprates: the 
magnitude of the PG decreases with increase of hole con- 
centration |,|,||. 

Figure 1(a) shows a theoretical I{V) tunneling char- 
acteristic in a superconductor-insulator-normal metal 
(SIN) junction (fig. 6 in ref. ||^). In the tunneling leghne, 
it is expected that the I{V) curve at high positive (low 
negative) bias, depending on the normal resistance of 
the junction, lies somewhat below (above) the normal- 
state curve. In conventional superconductors (SCs), the 
Blonder-Tinkham-Klapwijk (BTK) predictions are veri- 
fied by tunneling experiments [^. However, in cuprates, 
the BTK theory is violated. Figure 1(b) shows the SC- 
insulator-SC (SIS) I(V) curve measured in an under- 
doped Bi2212 single crystal with = 83 K (fig. 1 in ref. 
0). In fig. 1(b), one can see that the I{V) curve at high 
positive (low negative) bias passes not below (above) the 
straight dash line j|] but far above (below) the line. This 
fact cannot be explained by the d-wave symmetry of the 
order parameter. To our knowledge, this question has 
never been raised in the literature before. This finding is 
the main motivation of the present work. 

To the best of our knowledge, the soliton SC was for 
the first time considered in ref. ^ in order to explain the 
SC in organic quasi-one-dimensional (quasi-lD) conduc- 
tors. Later, Davydov JlO| , pl| applied the model of soliton 
SC to cuprates. The theory is based on the concept of 
bisolitons, or electron (or hole) pairs coupled in a singlet 
state due to local deformation of the lattice. 

Tunneling spectroscopy is an unique probe of SC state 
in that it can, in principle, reveal the quasiparticle 
(QP) excitation density of states (DOS) directly with 
high energy resolution. In this paper, we present tun- 
neling measurements performed on heavily underdoped 
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FIG. 1. (a) Theoretical I{V) tunneling characteristic in a 
SIN junction of a SC with the isotropic energy gap The 
dash line shows the normal-state curve, (b) Measured I{V) 
curve in a SIS junction of an underdoped Bi2212 with Tc = 
83 K [^. The dash line which is parallel to the I{V) curve at 
high bias is a guide to the eye. The arrows show the offset 
from the dash line. One can immediately notice the difference 
between the two plots pi. 



and slightly overdoped Bi2212 single crystals by using a 
break-junction technique. We find that m -plane tunnel- 
ing spectra below Tc are the combination of incoherent 
part from the PG and coherent QP peaks. There is a 
correlation between the magnitude of the PG and the 
magnitude of the SG in Bi2212. We find that the QP 
conductance peaks are caused by condensed topological 
solitons. Thus, we present an evidence for a quasi-lD 
topological-excitation liquid in Bi2212. 

The overdoped Bi2212 single crystals were grown us- 
ing a self-flux method as described elsewhere [|l^. The 
underdoped samples were obtained from the overdoped 
single crystals by annealing them in vacuum. The Tc 
value was determined by the four-contact method. The 
transition width is less than 1 K in the overdoped crys- 
tals, and a few degrees in the underdoped Bi2212. 

Experimental details of our break-junction setup can 
be found elsewhere p^ . In short, many break-junctions 
were prepared by gluing a sample with epoxy on a flexi- 
ble insulating substrate, and then were broken in the ab- 
plane by bending the substrate with a differential screw 
at low temperature in a He ambient. The electrical con- 
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FIG. 2. SIS dI/dV{V) and I{V) measured at 14 K within 
the same underdoped Bi2212 single crystal with Tc = 51 K. 
The dl/dViV) in both plots are normalized at -400 mV. The 
inset in the plot (b) shows the differences {dl / dV)a~{dI / dV)b 
and Ia,norm — Ib,norm- The inset in the plot (a) shows how la 
and lb are normalized: la is normalized at -400 mV, and 
is adjusted to be parallel at high bias to la (such procedure 
is equivalent to the normalization at ±cxd). 

tacts (typically with the resistance of a few ohms) are 
made by attaching gold wires to a crystal with silver 
paint. The I{V) and dI/dV{V) characteristics are deter- 
mined by the four-terminal method by using a standard 
lock-in modulation technique. At low (constant) temper- 
ature, in one junction, we usually obtain a few tunneling 
spectra by changing the distance between broken parts 
of a crystal, going back and forth etc., and, every time, 
the tunneling occurs most likely in different places. 

Figure 2(a) shows the SIS dI/dV{V) and !{¥) ob- 
tained in an underdoped Bi2212 single crystal, which look 
like usual spectra in Bi2212 |l^]. In fig. 2(a), the Joseph- 
son IcRn product is estimated to be 13.4 mV. The gap 
magnitude, Asc= 64 meV, is in good agreement with 
other tunneling measurements |jl^. The dI/dV{V) and 
I(V) shown in fig. 2(b) are obtained within the same 
underdoped single crystal as those in fig. 2(a). In fig. 
2(b), the gap having A — 130 meV is too large to be a 
SG. It is suggestive that the spectra in fig. 2(b) corre- 
spond to the PG. The SC in Bi2212 is weak in the heavily 
underdoped region |]l5| , |l6| . This may explain why it is 
possible to observe separately the PG in the heavily un- 
derdoped Bi2212 by taking into account that tunneling 
spectroscopy probes the local DOS. The absence of the 
Josephson current in the spectra shown in fig. 2(b) indi- 
cates that the humps in the conductance are incoherent. 
The differences between the spectra shown in figs 2(a) 




FIG. 3. SIS dI/dV{V) and I{V) measured at 14 K within 
the same overdoped Bi2212 single crystal with Tc = 88 K. The 
inset in the plot (a) shows SIN dl /dV{V) and I{V) measured 
at 9 K in an overdoped Bi2212 with Tc = 87.5 K (the same 
axis parameters as main plot). The dash lines in the plot (a) 
and in the inset, which are parallel at high bias to the I{V) 
curves, are guides to the eye. The inset in the plot (b) shows 
the normalized characteristics of a bound state of two solitons 
(a bisoliton) The gap is shown schematically in grey 

(the I{V) of the gap is not shown). The hight of conductance 
pea ks depends on the density of added (removed) electrons 



and 2(b) are presented in the inset of fig. 2(b), which 
correspond to a "pure SG" . Some parts of the dI/dV{V) 
in the inset of fig. 2(b) are slightly below zero because the 
spectra in figs 2(a) and 2(b) are not taken under the ex- 
act same conditions. The small humps in the dI/dV{V) 
shown in the inset of fig. 2(b) are discussed below. The 
dI/dV{V) and I{V) curves in the inset of fig. 2(b) re- 
semble the characteristics of a bound state of two solitons 
(a bi-soliton), shown in the inset of fig. 3(b). Thus, we 
find that the conductance in fig. 2(a) consists of the two 
contributions: the humps which correspond to the PG 
(which is a normal-state gap ||l|,||,||), and the QP peaks 
which are caused by condensed solitonlike excitations. In 
fig. 2(b), one can see that the PG is anisotropic. 

Figures 3(a) and 3(b) show the SIS dI/dV{V) and 
I{V) obtained within the same overdoped Bi2212 sin- 
gle crystal. In figs 3(a) and 3(b), the Josephson IcRn 
product is estimated to be 6 and 7.5 mV, respectively. 
The characteristics in fig. 3(a) look like usual spectra 
in Bi2212 The dI/dV{V) and I{V) curves in fig. 

3(b) resemble the spectra in the inset of fig. 2(b), and 
the bisoliton characteristics shown in the inset of fig. 
3(b). This means that the contribution from the PG 
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in the spectra shown in fig. 3(b) is small in compari- 
son with the contribution from the QP peaks, at least, 
at low bias. This is most likely due to the fact that, in 
slightly overdoped cuprates, the SC is the strongest, and 
the "strength" of the PG is weak [^,[6). At high bias, 
the contribution from the PG will be always predomi- 
nant, even, if the PG is weak. 

We also performed measurements in Ni-doped Bi2212 
single crystals (overdoped in oxygen): the measured data 
which are presented in an extended paper [ pTf are sim- 
ilar to the data in the overdoped Bi2212, shown in fig. 
3. From the temperature dependence of dI/dV{V), it is 
clear that the small humps which appear in dI/dV(V) 
at bias twice as large as the bias of the QP peaks, shown 
in fig. 3(b), relate to the QP peaks and not to the PG. 
The humps are also observed in the Ni-doped Bi2212. In 
the inset of fig. 2(b), a similar hump is present in the 
conductance at negative bias. These humps can be well 
understood in terms of a nanopteron soliton which is 
discussed below. 

Since the spectra measured in an underdoped Bi2212, 
shown in the inset of fig. 2(b), and the spectra measured 
in an overdoped Bi2212 (and Ni-doped Bi2212), shown 
in fig. 3(b), are similar, the data obtained in underdoped 
and overdoped Bi2212 are consistent with each other. 

In order to be sure that we observe not a SIS-junction 
effect but an intrinsic effect, we performed measurements 
in the overdoped Bi2212 crystals by SIN junctions. Pt- 
Ir wires sharpened mechanically are used as normal tips. 
The inset of fig. 3(a) shows the SIN dI/dV{V) and I[V) 
obtained in an overdoped Bi2212. In fig. 3(a), one can 
see that, basically, there is no difference between the I{V) 
characteristics measured in SIS and SIN junctions [see 
the dash lines in fig. 3(a) and in the inset of fig. 3(a)]. 

First, we give a description of a topological soliton. 
The topological soliton is an extremely stable nonlin- 
ear excitation which can be moving or entirely static 
|p8| , ^j20| . The solitons have particlelike properties, and, 
in solids, they occupy the intragap states in a CDW or 
SDW gap. The inset of fig. 3(b) shows the bisoliton 
characteristics. Since we consider a general solution, the 
gap shown schematically in the inset of fig. 3(b) is either 
a CDW or SDW gap. 

We now compare the measured data with theory. Soli- 
ton and bisoliton characteristics are described by hyper- 
bolic functions |p^ , p^ , pO| . Since the bisoliton conduc- 
tance peaks shown in the inset of fig. 3(b) look very 
similar to the conductance peaks not only of high- Tc SCs 
but also of low-Tc SCs (not the background) , we rely here 
exclusively on the I{V) characteristics which are concep- 
tually different for the two models: the BTK model for 
3D case and the model based on the concept of quasi-lD 
topological excitations. 

Figure 4(a) shows the measured I{V) curve from the 
inset of fig. 2(b). In fig. 4, for simplicity, we analyze 
the data only at positive bias. As shown in fig. 4(a), 
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FIG. 4. Measured /(V) curves and the f{V) fit (see text): 
(a) The measured data (diamonds) from tlie inset of fig. 2(b), 
the I{V) of the PG from fig. 2(b), and the f{V) fit. (b) The 
data from fig. 1(b), the f{V) fit, and their difference, (c) 
The I{V) from fig. 3(b), the f{V) fit, and their difference, 
(d) The SIN KV) from the inset of fig. 3(a), the fit (see text) 
and their difference. The Josephson currents are removed. 



the data from the inset of fig. 2(b) can be fitted very 
well by the hyperbolic function f{V) = Ax{tanh[{eV - 
2A)/eVo]+tanh[{eV + 2A)/eyo]), where e is the electron 
charge; V is the bias; A is the maximum SC energy gap, 
and A and Vq are the constants. In fig. 4(a), we also 
present the measured I{V) of the PG from fig. 2(b). We 
find that any I{V) characteristic obtained in Bi2212 can 
be resolved into the two components shown in fig. 4(a): 
from the quasi-lD topological excitations and from the 
PG. The "usual" I{V) and dI/dV{V) spectra in Bi2212 
show the presence of both components [see figs 2(a) and 
3(a)]. The absence [see fig. 2(b)] or weak contribution 
of one component [see fig. 3(b)] in spectra makes the 
appearance of the spectra " unusual" . 

Figure 4(b) shows the data from ref. 0, the f[V) fit, 
and their difference. Figure 4(c) depicts the two com- 
ponents in the I{V) curve from fig. 3(b). As shown 
in fig. 4(d), the contribution from the QP peaks in the 
SIN I{y) from the inset of fig. 3(a) seems to be weaker 
than that in SIS junctions. As seen in fig. 4, all plots 
are similar. To fit the SIN I{V), we use the same f{V) 
function by substituting 2 A for A [^^. In figs 4(b)- 
4(d), the amplitude. A, of the f{V) fit can be changed, 
this only affects the scale but not the shape of the dif- 
ferences which correspond to the PG. The I{V) curves 
in fig. 3(a) and in refs and jl^ can be resolved into 
the two components in the same manner. So, we con- 
clude that, in Bi2212, the I{y) characteristics of the QP 
peaks definitely disagree with the BTK theory, and are in 
good agreement with the theory of quasi- ID topological 
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The coherent QP peaks in dI/dV{V) curves can be fit- 
ted very weh by the derivative [/(V^)]' = Aix{{sech[{eV 
- 2A)/eVb])2 + isech[{eV + 2A)/eVo]f) 0. The fit is 
appUcable only to the QP peaks, and not to the humps 
at high bias, which correspond to the PG. 

In figs 2(a) and 3(a), the ratio of hump bias to QP-peak 
bias in the conductance remains more or less constant 
with doping, 2 (underdoped)-2.3 (overdoped), which is 
in good agreement with other tunneling data p4| . Since 
the humps correspond to the PG, the magnitudes of the 
PG and SG in Bi2212 correlate with each other. 

We now focus on the small humps which appear in 
dI/dV{V) outside the gap structure [see, e.g., fig. 3(b)]. 
They can be well understood in terms of a nanopteron 
soliton which is a bound state resulting from the nonlin- 
ear interaction between the soliton and the periodic wave 
(see fig. 6.27(c) in ref. @). Instead of flat asymptotics, 
the characteristics of a nanopteron have oscillations. The 
bisoliton can also exist in resonance with small ampli- 
tude linear waves |2^. The oscillations in asymptotics 
can be also seen in the I{V) shown in figs 4(a) and 4(c). 
What kind of periodic waves in Bi2212 can interact with 
the topological excitations and cause these oscillations? 
Phonons are the primary candidate on this role. On the 
other hand, since the magnetic resonance peak observed 
in inelastic neutron scattering spectra j2^ is caused by a 
magnon-like excitation (spin 1, charge 0), it is most likely 
that this spin excitation is in resonance with the bound 
state of topological solitons. 

It is well known that dI/dV{V) obtained in SIS junc- 
tions often have subgap structures [see, e.g., fig. 3(a)]. 
These subgap structures can be also understood in terms 
of a nanopteron soliton. They have most likely the same 
origin as the oscillations at high bias, i.e. due to the non- 
linear interaction between the bisoliton and the periodic 
wave (see fig. 4 in ref. ||2^ ). 

We now discuss the PG. In general, quasi-lD topologi- 
cal excitations reside on quasi-lD "objects". Recently, 
quasi- ID charge stripes separated by 2D AF domains 
have been discovered in some cuprates [^,^,^,^ . The 
presence of soliton excitations on stripes implies that the 
stripes are insulating, as predicted [p9[| , and not conduct- 
ing as widely believed. In other words, there is a CDW 
gap along stripes. Then, in order to pull out an elec- 
tron (a hole) from a stripe, it is necessary, first of all, 
to overcome the CDW gap on the stripe. Consequently, 
the tunneling PG is predominantly the CDW gap. The 
same conclusion can be independently obtained from the 
analysis of the data: The magnitude of PG, Aps= 130 
meV, shown in fig. 2(b), is too large to be explained by 
the development of local AF correlations, since the value 
of the superexchange energy, J ~ 120 meV [|3^, is not 
large enough to fit the data. Moreover, the magnitude 
of PG becomes larger than 130 meV at lower doping. 
Consequently, the tunneling PG is most likely a CDW 



gap. 

It is important to emphasize that the PG in transport 
measurements [ ]3l| , p^ is different from the tunneling PG. 
In transport measurements, the PG relates to a spin gap 
into local AF domains [^. As an example, if a molecular 
chain is embedded into a medium, a frictional force acts 
on the soliton [|ll| . 

We now turn to the model of SC in cuprates. The 
"fiat" (constant) asymptotics oi I{V) characteristics are 
the fingerprints of one-dimensionality. To the best of 
our knowledge, there are only two theoretical models of 
SC in cuprates which are based on the presence of one- 
dimensionality: the bisoliton model 1 10 1^ and the stripe 
model |3^. In the stripe model, charge stripes are as- 
sumed to be metallic [^, thus contrary to the experi- 
ment. By contrast, the bisoliton model is quantitatively 
in good agreement with the data: the bisoliton model 
predicts that by increasing the hole concentration the 
magnitude of the pairing gap decreases (see fig. 3 in ref. 
[Q). However, the bisoliton model is a theory of soliton 
pairing, but the mechanism of the establishment of phase 
coherence is not considered. Experimentally, spin fluctu- 
ations mediate the phase coherence in cuprates [Q . The 
bisolitons are formed due to electron-phonon interactions 
which are moderately strong and nonlinear [HO 1 1 1 . 



Lastly, it is worth noting that the analysis of many 
data measured in cuprates shows that the data can be 
naturally understood in the framework of the quasi-lD 
topological-excitation-liquid scenario [|l^ . 

In summary, tunneling measurements have been car- 
ried out on underdoped and overdoped Bi2212 single 
crystals. Tunneling spectra below Tc are the combina- 
tion of incoherent part from the pseudogap and coherent 
quasiparticle peaks. There is a correlation between the 
maximum magnitude of the pseudogap and the distance 
between the quasiparticle peaks. The I{V) characteris- 
tics of the quasiparticle peaks are in good agreement with 
the theory of quasi- ID topological excitations. The soli- 
tons and bisolitons reside most likely on charge stripes. 
It seems that magnon-like excitations which cause the 
appearance of the magnetic resonance peak in inelas- 
tic neutron scattering spectra are in resonance with the 
bisolitons. The bisoliton model of superconductivity in 
cuprates |l0| , |rT| is correct in the description of pairing 
characteristics, but it lacks the mechanism of the estab- 
lishment of phase coherence. 
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